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Abstract 

We explore strong gravity effects of the geodesic motion in the spacetime of rotating black 
' holes in general relativity and braneworld gravity. We focus on the description of the motion in 
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' the innermost stable circular orbits and beyond them as well as at the characteristic stable orbits, 
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terms of three fundamental frequencies: The orbital frequency, the radial and vertical epicyclic 
frequencies. For a Kerr black hole, we perform a detailed numerical analysis of these frequencies at 



at which the radial epicyclic frequency attains its highest value. We find that the values of the 
epicyclic frequencies for a class of stable orbits exhibit good qualitative agreement with the observed 



\ frequencies of the twin peaks quasiperiodic oscillations (QPOs) in some black hole binaries. We 

also find that at the characteristic stable circular orbits, where the radial (or the vertical) epicyclic 
frequency has maxima, the vertical and radial epicyclic frequencies exhibit an approximate 2 : 1 
ratio even in the case of near-extreme rotation of the black hole. Next, we perform a similar 
analysis of the fundamental frequencies for a rotating braneworld black hole and argue that the 
existence of such a black hole with a negative tidal charge, whose angular momentum exceeds the 
Kerr bound in general relativity, does not confront with the observations of high frequency QPOs. 
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I. INTRODUCTION 



Black holes are one of the most exciting objects of study in modern theoretical physics 
and astrophysics. They occupy a central place in all theories of gravity formulated in various 
spacetime dimensions. 



A. Black holes in general relativity and braneworld gravity 

In four dimensions, general relativity (GR) admits a family of stationary black hole 
solutions which turns out to be crucial for understanding its nature and the occurrence of 
spacetime singularities py]. These solutions possess a number of striking properties such as 



uniouene. and ^tabiUty, Udden symmetnes and integ^abUity of geode.c= etc. M (- a..o 
[5|, |6[). Altogether, these properties pave the way for astrophysical implications of GR in 
the regime of strong gravity, stimulating the search for black holes in the physical universe. 

An astrophysical black hole is thought to be described by an exact solution of GR, 
discovered by Kerr W. This solution is uniquely characterized by two parameters: The 
mass and the angular momentum. In ^, it was first shown that the Hamilton- Jacobi 
equation for geodesies in the Kerr metric admits a complete separation of variables that is, 
the geodesic motion is completely integrable in this spacetime. In subsequent developments, 
various aspects of the geodesic motion in the Kerr metric have been extensively explored 
by many authors (for instance, see 5[ and references therein). Among these results, the 
calculations of the observable orbits of test particles are of particular interest. The authors 
of js] were the first to give a complete description of these orbits for the circular motion in 
the equatorial plane of the Kerr black hole. In particular, they found that for an extreme 
rotating Kerr black hole, the binding energy of the innermost stable circular orbit (ISCO) 
(the maximum amount of energy released by a particle approaching this orbit) can attain 
nearly 42% of the particle rest-mass energy, whereas for a Schwarzschild black hole it is 
about 6%. The high efficiency of this process led to the idea of invoking an accretion disk 
around a black hole to explain enormous energy output from both X-ray binaries and active 



galactic nuclei (for details, see {q]). 

In higher dimensions, there exist a number of important black hole solutions with un- 
expected physical characteristics, such as dynamical instability, different horizon topologies 



and different rotation dynamics. Among them, the most interesting blaclc hole solution is 
given by the Myers-Perry metric jl^, which is supposed to be a generalization of the Kerr 
solution to all number of spacetime dimensions. However, this solution is not unique, unlike 
the Kerr solution. In five-dimensional vacuum gravity there exists a rotating black ring 
solution [11], which may have the same mass and angular momentum as the Myers- Perry 
solution. 

On the other hand, the higher-dimensional black holes may have observable consequences 
as well. We recall that some phenomenological aspects of string/M theory lead to braneworld 
scenarios, which imply that our physical universe is geometrically described by a (3-1-1)- 
dimensional hypersurface (a "3-brane") embedded in higher dimensions 12|, ll3|]. A com- 
plete description of black holes in such scenarios is a thorny issue and many aspects of it 
still remain open. Remarkably, there are several simplifying approaches which provide an 
intriguing description. For instance, if the size of the black hole is much smaller than the 
scale of the extra dimensions, the black hole would behave as a higher-dimensional object. 
To good enough approximation, one can argue that these black holes are described by Myers- 
Perry type solutions of the Einstein field equations in higher dimensions. Furthermore, the 
fact that in the braneworld scenarios the scale of quantum gravity becomes as lower as 
TeV-energy scales indicates that the small-size black holes would have been detected at high 
energy experiments [l4, 15] . Of course, there always exists the question of gravitational 



stability of the Myers-Perry type black holes on the brane, whose final status is still not 



16|. 



clear. Some interesting aspects of this issue were investigated in 

As for the large black holes, for which the size of the horizon is much larger than the length 
scale of the extra dimensions, they would also carry the imprints of the extra dimensions 
though "effectively" they must look like four-dimensional objects. In one approach employed 



m 



ITj . it was suggested that such a black hole in the Randall and Sundrum (RS) braneworld 
scenario could be described by the ordinary Schwarzschild solution on the brane, which in 
five dimensions would look like a black string. The main drawback here is that the black 
string solution suffers from curvature singularities, propagating along the extra dimension. 



This solution also suffers from classical Gregory-Lafiamme instabilities 18|. In another 



approach, such black holes were described by postulating the metric form and solving the 



effective gravitational equations on the RS brane 19|, |20|. This idea is also reinforced by 



the fact that in the low energy limit, the RS braneworld scenario supports the properties of 



four- dimensional GR on the brane [2l|, |2^. It is interesting that in this approach one can 
analytically describe, to some extent, the braneworld black holes by a Reissner-Nordstrom 
;ype solution in the static case and by a Kerr-Newman type solution in the rotating case 



20| . The tidal charge appearing in both solutions and being the imprint of the extra fifth 



dimension supersedes the square of the electric charge in the usual Reissner-Nordstrom and 
Kerr-Newman solutions. Furthermore, the tidal charge in these solutions may have both 
positive and negative values. On a par with the Kerr black hole, the rotating braneworld 
black hole can be considered as the useful and the simplest model to explore the physical 
effects of the tidal charge, though the gravitational stability of the latter remains an open 
problem. Some astrophysical implications of the braneworld black holes have recently been 



discussed in 



23|-|26|. 



B. High-frequency quasiperiodic oscillations 

Nowadays, astronomical observations provide compelling evidence for the existence of 



black holes in the universe. To focus only on black holes in X-ray binaries, today 23 b 



ack 

3- 



holes are known and 32 X-ray binaries are thought to contain black hole candidates 
Observations of the X-ray binaries with both black holes and black hole candidates have 
revealed finite-width peaks in the X-ray spectra, which can be thought of as signatures of 
high-frequency (> 40 Hz) quasiperiodic oscillations (QPOs) of the black hole accretion disk. 
Furthermore, in a number of sources, two peaks of QPOs have been detected. The associated 
frequencies of these peaks are (300, 450 Hz) for X-ray binary GRO J1655-40; (113, 168 Hz) 
for X-ray binary GRS 1915 + 105; (184, 276 Hz) for X-ray binary XTE J1550-564 and (165, 
241 Hz) for X-ray binary H 1743-322, harboring a black hole candidate. The twin frequencies 
are in a 3:2 ratio. The source GRS 1915 + 105 also exhibits a pair (164, 328 Hz) which is 



in a 2:1 ratio (see 



28 



for further details). 

Observations of such frequencies are certainly of fundamental importance as they may 
provide crucial probes of strong gravity near the black holes. Moreover, when combined with 



the measurements of the angular momenta of black holes [29 



they raise once again the 



old compelling question: What is the bona fide geometry of the spacetime around the black 



holes'! For instance, as argued in 30| , the lower bound on the angular momentum of a black 



hole in the X-ray binary GRS 1915 + 105 turns out to be very close to the (upper) Kerr 



bound in GR. This raises the question: Does the Kerr solution of GR describe the observed 
black holes with the high value of the angular momentuml Today it remains unknown how 
future observational data will deal with these questions. However, this perspective in itself 
is very challenging for further study of observable effects of black holes in both GR and 
beyond it, including the cases with the violation of the Kerr bound. Some effects of rotating 
)lack holes in string theory, for which the Kerr bound is breached, were recently studied in 



3l|. 



It is also interesting that in a number of theoretical models explaining the origin of 
high frequency QPOs, general relativistic epicyclic frequencies of the geodesic motion have 
3een considered as " building blocks" , though the precise mechanism still remains unknown 
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32 



34| . The successive theory of the epicyclic motion around the Kerr black holes 



in general relativity (with an electric charge or with an external uniform magnetic field) 
was developed in 35| (see also companion papers 36|-|38|). The authors, for the first time, 
obtained closed analytical expressions for the epicyclic frequencies of radial and vertical 
oscillations and put forward an idea of nonlinear resonances, which may occur when these 



frequencies are in a rational relation 35| . Afterwards, ana^ 



of radial and vertical oscillations were also appeared in 



39 



ytic expressions for the frequency 



40| in the context of the analysis 



of trapped oscillations of an accretion disk around the Kerr black hole. The theory of 



epicyclic motion and the idea of nonlinear resonances have been utilized in 33 



reproducing in some cases the plots of 



34j . even 



35j for the positions of nonlinear resonances, to 



explain the observed frequency ratios of high frequency QPOs in black hole binaries. 



C. Motivation and the basic idea 

In light of all these developments, it becomes of interest to return once more to the 
exploration of the observable effects of a strong gravity regime encoded in the geodesic motion 
around rotating black holes. Our motivation for this arises from the fact that the epicyclic 
frequencies of the geodesic motion, as mentioned above, lie at the root of theoretical models 
for the observed high frequency QPOs in black hole binaries. Furthermore, the continuing 
success with the measurements of the angular momentum of black holes motivates us to 
explore the effects of the geodesic motion not only for rotating black holes in GR, where 
the Kerr bound on the angular momentum holds, but also for rotating black holes in a 
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braneworld, for which the Kerr bound may be breached. The basic idea of the present paper 
is as follows: 

(i) To perform a detailed numerical analysis of the epicyclic frequencies at the stable 
circular orbits of interest in the Kerr field, using the general analytical expressions found 



in earlier papers 



35 



36|. Thus, "monitoring" the behavior of these frequencies at the 
stable orbits, especially at those where the radial (or the vertical) epicyclic frequency attains 
its highest value, we intend to clarify the question: Whether the simple model of pure 
general relativistic geodesic motion can provide observationally viable values for the epicyclic 
frequencies and their ratios'? 

(ii) To extend the numerical analysis of the epicyclic frequencies to the case of a rotating 
braneworld black hole as well. The fact that the associated spacetime metric has the Kerr- 
Newman type form with the tidal charge (instead of the square of the electric charge) and for 
the negative tidal charge, the angular momentum of the black hole exceeds the Kerr bound 
makes this model very instructive. It is also important to note that the circular motion in 
this model may occur in stable orbits, unlike the case of the ordinary black holes in higher 
dimensions, for instance, the case of Myers-Perry black holes [41]. Exploring the behavior of 
the epicyclic frequencies at particular orbits of interest and comparing the results with those 
for the Kerr black hole we want to address the issues: How the presence of the tidal (positive 
or negative) charge influences the values of the epicyclic frequencies? Does the potential 
existence of over-rotating braneworld black holes confront with modern observations of high 
frequency QPOs in black hole binaries'? 

The organization of the paper is as follows: In Sec. II we begin by a brief review of the 
salient features of the geodesic motion in the Kerr metric and give analytic expressions for 
the epicyclic frequencies of radial and vertical oscillations. Next, we present the results of 
a detailed numerical analysis of these frequencies at the stable circular orbits of interest. 
In Sec. Ill we briefly discuss the salient properties of a rotating black hole, localized on a 
3-brane in the Randall-Sundrum scenario as well as those of the geodesic motion in the field 
of this black hole. Here we present analytic expressions for the epicyclic frequencies of radial 
and vertical oscillations and perform their full numerical analysis, focusing on the stable 
orbits, where the radial epicyclic frequency reaches its maximum value. We also discuss the 
observational signature of the negative tidal charge carried by the braneworld black hole 
and argue that the idea that such a black hole can potentially exist in the universe does not 



confront with the observations of high frequency QPOs in black holes binaries. In Sec. IV 
we conclude with a brief discussion of our new results. 



II. GEODESIC MOTION AND EPICYCLIC FREQUENCIES IN THE 
SPACETIME OF KERR BLACK HOLES 

We begin with the Kerr spacetime which in the Boyer-Lindquist coordinates is given by 
the metric 



9 ^ f dr^ ,^o\ sin^6' 



ds [dt- a sin^ 9 d(j)] + [ — + dO'' \ + adt - (r'' + a^) d(f) , (1) 

where 

A = r^ + a^ - 2Mr , S = + cos^ 9 , (2) 

and M is the mass, a is the rotation parameter or the angular momentum per unit mass, 
a = J/M. The most striking feature of this spacetime is that it contains a rotating black 
hole. This becomes evident by the existence of the event horizon which is governed by the 
equation A = 0. The largest root of this equation given by 



r+ = M + VM2 - a2 (3) 

determines the radius of the horizon. The event horizon, as follows from this expression, 
exists provided that a < M . Thus, a rotating black hole in general relativity must possess 
an angular momentum not exceeding its mass. 

Another remarkable feature of the Kerr spacetime is that it admits the complete separa- 
tion of variables in the Hamilton- Jacobi equation for geodesies. The global time-translational 
and rotational symmetries of this spacetime, along with its hidden symmetries, make the 
Hamilton- Jacobi equation completely integrable ji] (see also {sl]). For our purposes in the 
following, we will focus only on the motion of test particles in the equatorial plane 6 = it/2. 



From the symmetry considerations, it 
(cyclic) orbits. Following works of 35 



bl 



ows that such a motion will occur in circular 
one can simply describe the circular motion 
and its perturbations, i.e. off-equatorial (epicyclic) motion, using the method of successive 
approximations in the geodesic equation 
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where the parameter s denotes the proper time along the geodesies and V^f^ are the Christof- 
fel symbols of the spacetime metric. Clearly, for the circular motion in the equatorial plane 
r = Vq, 9 = 7r/2, one can introduce the position vector = {t{s) ,ro ,vr/2 , f2ot(s)}, where 
Qq is the orbital frequency. Next, using the Christoffel symbols for the Kerr metric it is easy 
to show that for the circular motion, the /i = 0,2,3 components of equation (jlj) become 
trivial, whereas the remaining component with fi = 1 yields 

where fig = M^/^/r^/^ is the usual Kepler frequency, the upper sign refers to direct orbits 
(the motion of the particle and the rotation of the black hole occur in the same direction) 
and the lower sign refers to retrograde orbits (the motion of the particle is opposite to 
the rotation of the black hole). Using equation (jS]) in the normalization condition for the 
four- velocity, g^yU^u^ = — 1 , we obtain the expression for energy of the particle 



E _ - 2Mr ± a VMr 

m r {r^ -3Mr ±2aVM^y^^ ' 

where E = muQ. This expression was first found and studied in jsj. The vanishing denom- 
inator of this expression, determines the radius Vph of the limiting photon orbit. That is, 
the circular motion exists only in the region r > r^/j. It is not difficult to show that for 
an extreme Kerr black hole, a = M, we have = M for direct orbits and r^/j = 4M for 
retrograde orbits, while for a = 0, we find that Vph = 3M. 

To describe the epicyclic motion, we introduce a deviation vector 

e{s)=x'^{s)-z'^{s), (7) 

and expand equation (jlj) in powers of C,^{s) about the circular orbits. Focusing on the 
linear approximation, it is straightforward to show that the epicyclic motion amounts to 
two decoupled oscillations in the radial and vertical directions (the details can be found in 



35l . |36| and in recent works 25|, |42|). For the frequency of radial oscillations, we have 



6M 3a 



2 



Qi = Q'o {l- — - — ±8aQA, (8) 



whereas the frequency of vertical oscillations is given by 



nl = nl [1 + — T 4a^], . (9) 



The stability of the circular motion against small oscillations is determined by the conditions 
Q^. > and Qg > 0. From the condition for the radial stability, one can infer the radii of 
the innermost stable circular orbits (ISCOs). As it was first shown in [8j, for a nonrotating 
black hole, a = 0, we have r^s = 6M, while for the extreme rotating M, we find 

that = M for the direct orbit and r^s = 9M for the retrograde one. Meanwhile, it is 
easy to verify that the expression in (|9]) is always nonnegative in the region of existence and 
radial stability of the circular motion. That is, the motion is stable with respect to small 
oscillations in the vertical direction. 

To summarize, the description of the circular motion and its small perturbations in the 
Kerr field inevitably results in three different fundamental frequencies, Qq, Qr and Qg , unlike 
the case of Newtonian gravity, where all the three frequencies turn out to be the same as 
the Kepler frequency, Qq = Qr = = ■ (We note that for the Schwarzschild black hole 
two of these different frequencies are different. Namely, Qr and Qg = Qq = Qg). It is the 
regime of strong gravity near the Kerr black hole, combined with its rotational dynamics, 
that distinguishes between the three fundamental frequencies. Therefore, the values of 
these frequencies at the radii of physical interest are of great importance for astrophysical 
implications of the Kerr black hole. 

In the following, we use the coordinate frequency ly in physical units, instead of the 
angular frequency Q = 27rz/, as well as the characteristic frequency scale ui = c' /2t{GM ~ 
3.2 ■ IQ'^ {Mq/M) Hz . We recall that here c is the speed of light, G is the gravitational 
constant and Mq is the mass of the Sun. With this in mind and using expression ([5]), we 
find that the orbital frequency around an extreme rotating black hole (a = M) is 

uo = \u,^ 1.6 ■ 10^ (^) Hz (10) 

for the direct ISCO whereas, it is given by 

''°-4^'-^-2-i'''(§)"^ 

for the retrograde ISCO. We have performed a detailed numerical analysis of the orbital 
and vertical frequencies at ISCOs around a Kerr black hole with mass M = IOMq. The 
results of calculations are given in Table I. We note that as the rotation parameter of the 
black hole grows, the radius of the direct ISCO moves towards the event horizon and the 
associated orbital frequency increases, approaching its maximum value in (|T0|) . The vertical 

9 



epicyclic frequency attains its maximum value and then decreases to zero, i.e. z/g = for 
a = M, r = M. The ratio of the frequencies, z/o/z/g, essentially differs from unity only 
for the fast enough rotation of the black hole. For the retrograde motion both frequencies 
decrease with the growth of the rotation parameter (see also Eq. flTT]) ). whereas their ratio 
remains about unity. 

It is also of interest to calculate all the three frequencies at direct stable orbits beyond 
ISCO (r > rjsco) around an extreme rotating black hole, a = M. The results of numerical 
calculations are presented in Table II. We observe that at particular orbits near the black 
hole, the predicted values of the epicyclic frequencies (bolded in the table) are close to 
the corresponding frequencies of twin peaks QPOs, which are in a 3:2 ratio and have been 
detected in some black hole binaries: For instance, to (184, 276 Hz) for X-ray binary XTE 
J1550-564; to (165, 241 Hz) for X-ray binary H1743-322; to (113, 168 Hz) for X-ray binary 
GRS 1915 + 105. 



A. The highest epicyclic frequencies 

It is interesting that there exist the highest epicyclic frequencies of small oscillations 
around circular orbits in the Kerr field. This fact was first noted in 39], for the frequency of 
radial oscillations given by expression ([8]). Evaluating the first derivative of this expression 
with respect to r, we obtain the equation 

(8M -r) + a^ {5r^ - 4Mr) ± 2ay/Mr [a^ + r {M - 6r)] = . (12) 

This equation determines the radii for both direct and retrograde orbits, at which the radial 
epicyclic frequency attains its highest value. For the Schwarzschild = 0, it follows 

that r^ax = 8M and the associated frequency i^rimax) — 707.1 (Mq/M) Hz . Meanwhile, for 
a 7^ 0, equation fll2p can be solved only numerically. In particular, for a = M and for the 
direct orbit, we find that Vmax — 2.4M and 

iyr^max)^2A53(^)Hz. (13) 



M 

Similarly, for the retrograde orbits at a = M, we have r^ax — 11. 8M and 

'Mrr 



lyrimax) ^ 422.6 (^^) Hz . (14) 



10 



In Figure 1 we plot the dependence of the radial epicyclic frequency on the radii of 
circular orbits around the Kerr black hole for different values of the rotation parameter and 
for M = IOMq. The full numerical analysis of equation (IT^ and the associated values of 
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FIG. 1. Radial epicyclic frequencies with three values of the rotation parameter a/M = , 0.5 and 
0.99. (Left: For direct orbits. Right: For retrograde orbits). 

the radial, vertical and orbital frequencies along with their corresponding ratios are given in 
Tables III and IV. 

Comparing Figure 1 and Table III, we see that for the direct orbits, with increasing 
rotation parameter of the black hole, the maxima of the radial epicyclic frequency shifts 
towards the event horizon and in the limiting case a = M, the frequency attains its highest 
value in the near-horizon region. The accompanying vertical and orbital frequencies at the 
same radii also increase to their highest values for a = M. It is also interesting to note that 
the characteristic ratios : Vrimax) =2:1, z/q : Vrimax) =2:1 and Uq : uq = 1 : 1 remain 
almost unchanged up to large enough values of the rotation parameter. However, for a — )■ M 
we have the approximate ratios ug : Ur = 9 : 5, Uq : Ur = 5 : 2 and Uq : ug = 3 : 2 . Thus, we 
conclude that at the characteristic stable circular orbits, where the radial epicyclic frequency 
attains its highest value, the ratio ug : Uj. = 2 : 1 remains nearly the same even for a — )■ M. 
Remarkably, this fact is in good enough agreement with the observed twin QPOs frequencies 
in the X-ray spectrum of some black hole binaries. For instance, for a ^ (0.8-0.9)M, the 
detected pair (164, 328 Hz) in the source GRS 1915 + 105 falls in the expected ranges of 
the radial z/,. and vertical ug epicyclic frequencies given in Table IIIII Meanwhile, Figure 1 
and Table IV for the retrograde orbits show that the highest value of the radial frequency 
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at a = as well as the associated values of the vertical and orbital frequencies decrease with 
the growth of the rotation parameter and attain their characteristic values in the limiting 
case a = M. It is also interesting to note that the frequencies exhibit, to a good enough 
accuracy, the ratios : I'rimax) =2:1, Uq : Urimax) =2:1 and Uq : ug = 1 : 1 . 

We turn now to the expression of the vertical epicyclic frequency given in For direct 
orbits and for sufficiently large values of the rotation parameter, this frequency attains its 
highest value as well. The radii of characteristic orbits, pertaining to the maxima of the 
vertical frequency, are given by the equation 



+ (5r - 2M)1 + 2aVMT (a^ - 3r^) = . (15) 



r 

Solving this equation numerically for a = M, we find that r(^max) — 1-86M and the highest 
vertical frequency is 

i^eimax)^ ^875 (^^'^ Hz, (16) 
whereas, the radial frequency at this radius has the value 

Ur- 223Q(^^^Hz. (17) 

It is easy to see that the approximate ratio I'e^max) : z^r = 2 : 1 holds in this case as well. 

In Figure 2 we plot the vertical epicyclic frequency as a function of the radius of direct 
orbits, for given values of the rotation parameter and for M = IOMq [Left). Here we also 
plot the positions of ISCOs and ve{max) as the functions of the rotation parameter [Right). 
We see that the vertical frequency reaches its highest value in the region of physical interest, 

^{max) > f I Seo- 



ul. GEODESIC MOTION AND EPICYCLIC FREQUENCIES AROUND 
ROTATING BRANEWORLD BLACK HOLES 

The measurements of the angular momenta of black holes harbored in a number of X-ray 
binaries reveal high enough values, approaching in some cases the extreme Kerr limit. For 
instance, the lower hound on the rotation parameter of the black hole in the X-ray binary 
GRS 1915 + 105 turns out Jx) be very close to the {upper) Kerr bound for rotating black 
holes in general relativity [30]. Although similar results are somewhat model-dependent, 
nevertheless they reinforce the belief that a black hole that possesses an angular momentum 
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FIG. 2. Direct orbits: Vertical epicyclic frequencies with the rotation parameter a/M = 0,0.5 
and 0.99 (Left). Positions of ISCOs and ve{max) the functions of the rotation parameter (Right). 

greater than the Kerr bound could exist. This in turn stimulates the study of the observable 
effects of black holes with over-Kerr bound angular momenta ("over-rotating" black holes). 
In this section, we explore this possibility in an instructive model of geodesic motion around 
a rotating black hole localized on a 3-brane in the Randall- Sundrum scenario 20|. Such a 
black hole carries a tidal charge, as an imprint of the extra fifth dimension, and breaches 
the Kerr bound for the negative value of the tidal charge. The associated spacetime metric 
in Boyer-Lindquist coordinates is given by 



ds"^ = — — (dt — a sin^ 9 1 



+ 



sin^ 6 



adt 



where 



A = + - 2Mr + /3 , 



S = + cos^ 9 . 



(18) 



(19) 



Here M and a are the mass and rotation parameters, respectively and the parameter /3 is 
the tidal charge. This metric is the reminiscent of the usual Kerr-Newman metric in which 
the electric charge of the black hole is superseded by the tidal charge. Solving the equation 
A = 0, we find that the radius of the event horizon is given by 



r+ = M+ JA'P -a'^-13, 



which implies the relation 



(20) 



(21) 



13 



where the equahty holds in the extreme case. We see that for the positive tidal charge, the 
angular momentum of the braneworld black hole does not exceeds its mass, whereas for the 
negative tidal charge it is greater than the mass. The detailed analysis of the properties of 



metric f[T5]) can be found in 20l. |25|. 



The circular motion of test particles in metric (fT8|) occurs in the equatorial plane 6 = 7t/2, 
whose orbital frequency Qo has the same form as given in ([5]), where the Kepler frequency 
is now given by ils = {Mr — /r^ . As for the energy of the test particles, it can be 
found either by substituting in equation the nonvanishing components of the Christoffel 
symbols and repeating all the steps leading to equation (E]) or simply replacing the square 
of the electric charge by the tidal charge, in the corresponding expression for the usual 



Kerr-Newman metric (see, for instance 



351]). As a result, we have 



E 
m 



r2 - 2Mr + /3 ± a ^Mr - P 

172 



3Mr + 2;^ ± 2 a y/Mr - (3 



(22) 



The vanishing denominator of this expression gives the radius of the limiting photon 
orbit. For /3 = M^, the limiting radius is the same as for an extreme Reissner- Nordstrom 
black hole, rph = M (r+ = M). Meanwhile, for the negative tidal charge /3 = — 
(a = \/2M , r+ = M), the limiting radii are rph = M for the the direct motion and rph = 
4.82M for the retrograde one |2o|. 

As for the epicyclic motion, it is frequency in the radial direction is given by 



"0 



Mr- (3 

while, for the frequency in the vertical direction we have 

2Mr -/?■ 



(23) 



1 + 1 + 



Mr-/3 



2Mr - 13' 
' Mr- 13 



(24) 



We note that these expressions were also given in 25| and they can be obtained from the 



general expressions for the Kerr-Newman field, earlier found in [35|, in which one needs to 
change the square of the electric charge to the tidal charge, keeping in mind that the latter 
may have both positive and negative values. For the vanishing tidal charge, /3 = 0, these 
expressions go over into those given in (jH]) and ([9]). The boundaries of the ISCOs in the 
radial direction are governed by the equation fi^ = . The analysis of this equation shows 
that for /3 = M^, the limiting radius is r^s = 4M just as for an extreme Reissner- Nordstrom 
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black hole, whereas for an over-rotating black hole with [5 = — (a = \/2M)we 



have 



rms = M for the direct motion and r^s — 11.25M for the retrograde one (see [20|, l25|] for 
details). Meanwhile, the analysis of expression shows that the vertical motion is stable 
against linear perturbations in this direction. 

It turns out that the radial epicyclic frequency in fl23|) has maxima at some characteristic 
orbits, just as expression ([S]) for the Kerr field. In what follows, we will focus on this case. 
Assuming first that the black hole has a small positive tidal charge, we compute all three 
fundamental frequencies and their corresponding ratios at direct orbits, for which the radial 
frequency in (123!) attains its maximum value. The numerical results are summarized in 
Table |Vl 

Comparing these results with those given in Table IIIIl we see that the observed pair of 
frequencies (164, 328 Hz) in the source GRS 1915 + 105 falls in the range of the radial u,. and 
vertical frequencies that corresponds to less values of the rotation parameter, a ~ (0.7- 
0.8)M. On the other hand, recent observations give the lower bound a > 0.98M on the 

n 

rotation parameter of the black hole in GRS 1915 + 105 [30[. Nevertheless, this discrepancy 
still might not be enough on its own to judge about the observational appearance (or non- 
appearance) of the positive tidal charge as the reported constraint on the bound of the 
rotation parameter is highly mo del- dependent. In some other cases, this constraint is found 



to be much below (see, for instance. 
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44j|). The reason for this is that the measurements 
of the rotation parameter of the black hole are obtained by modeling the thermal X-ray 
continuum of its accretion disk spectra to infer the inner radius of disk with subsequently 
identifying it with that of ISCO, which depends only on the mass and rotation parameter of 
the black hole. However, the accretion disk spectra are generally complex being accompanied 
by a radiative tail of non-thermal nature as well as by the imprints of relativistic effects. 
Depending on the model under consideration, these may have a large impact on the estimate 
of the rotation parameter. 

Meanwhile, there exists a fundamental constraint on the values of the positive tidal 
charge and the rotation parameter, governed by equation (1201) . i.e. by the very existence 
of the black hole. It follows from this equation that with the growth of the positive tidal 
charge, the rotation parameter must decrease and for the large enough charge, /3 — )■ , the 
rotation parameter tends to zero. Furthermore, performing the above numerical analysis of 
the epicyclic frequencies for the tidal charge /3 = 0.99M^ , which implies the highest rotation 
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parameter a = O.IM , we find the pair of frequencies {vrii-'e )= (118, 249 Hz). This pair 
lies far enough away from the detected pair of frequencies (164, 328 Hz) in the source GRS 
1915 + 105. On the basis of all that has been said above, we can conclude that at least, the 
large enough positive tidal charge is not supported by the observations of QPOs in black hole 
systems. 

Next, we suppose that the black hole possesses the negative tidal charge (3 = — and 
again calculate all three frequencies at the characteristic direct orbits, at which the radial 
epicyclic frequency attains its maxima. The results are given in Table IVI[ Comparing these 
results with those given in Table IIIH it is easy to see that the negative tidal charge has an 
enlarging effect on the radii Vmax, whereas with increasing the rotation parameter of the 
black hole, the radii again move towards the event horizon, approaching the limiting value 
for a = y/2M. We also see the appearance of the approximate 2 : 1 ratio between the 
vertical and radial epicyclic frequencies. It is important to note that in the over-rotating 
case, a ~ 1.27M, the values of the radial Ur and vertical z/g frequencies and their ratio are 
in good enough agreement with the detected pair of frequencies (164, 328 Hz) in the source 
GRS 1915 + 105. 

Figure [3] displays the positions of the maxima of the radial epicyclic frequencies for an 
extreme Kerr black hole and for an over-rotating braneworld black hole with the negative 
tidal charge /3 = — M^. We see that the epicyclic frequencies in the field of these black holes 
are observationally almost indistinguishable. We note that it is the value of the angular 
momentum (for instance, reported from independent observations) that would play a crucial 
role to distinguish between these two black holes and the over-Kerr bound value of the 
angular momentum would be in favor of the potential existence of the braneworld black hole 
in the universe. 
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FIG. 3. The radial epicyclic frequencies for the direct orbits. The soUd Une corresponds to the 
extreme Kerr black hole, /3 = 0, a = M and the dotted line refers to the braneworld black hole, 
/3 = -M2, a = ^/2M. 

IV. CONCLUSION 

In general relativity, the circular motion of test particles around rotating black holes can 
be described in terms of three fundamental frequencies: The orbital (cyclic) frequency, the 
epicychc frequency of radial oscillations and the epicyclic frequency of vertical oscillations 
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38| . All three of these frequencies are different, encapsulating somewhat the effects of 



strong gravity regime near the black holes. Remarkably, the epicyclic frequencies of radial 
and vertical oscillations can be related to the frequencies of twin peaks high-frequency QPOs, 
discovered in a number of black hole binaries. The purpose of this paper was to explore 
these issues in more detail from the general relativistic point of view. We have considered 
two instructive models of the geodesic motion; around a Kerr black hole in general relativity 
and around a rotating black hole in braneworld gravity. In both cases, we have performed 
a detailed numerical analysis of the fundamental frequencies at the stable orbits of interest 
and found the examples of observationally viable cases of the epicyclic frequencies. 

As the main result for the Kerr black hole, we have found that for a class of stable 
orbits, where the radial (or the vertical) epicyclic frequency attains its maxima, the vertical 
and radial epicyclic frequencies are in an approximate 2 : 1 ratio, occurring even for near- 
extreme rotation of the black hole. We have also found that for the rotation parameter in the 



17 



range a ~ (0.8-0.9)M, the values of the epicychc frequencies underlying this ratio exhibit 
good enough qualitative agreement with the detected twin peaks frequencies of QPOs in 
the source GRS 1915 + 105 which, to our best knowledge, appears to be a unique case so 
far. Furthermore, we have found that at particular class of orbits near the black hole, the 
predicted values of the epicyclic frequencies (bolded in Table II) are close to the frequencies 
of twin peaks QPOs, which are in a 3:2 ratio and have been detected in some black hole 
binaries: For instance, to (184, 276 Hz) for X-ray binary XTE J1550-564; to (165, 241 Hz) for 
X-ray binary H1743-322; to (113, 168 Hz) for X-ray binary GRS 1915 + 105. These findings 
are substantially new. As we noted in the introduction, the idea of nonlinear resonances that 
may happen when the epicyclic frequencies are in a rational relation 35(] has been utilized 



in the literature [33|, 



34| to develop the models for the observed high frequency QPOs. The 



results of our paper show that the simple model of pure geodesic motion can also provide, to 
some extent, the observationally viable values for the epicyclic frequencies, without invoking 
nonlinear resonances between them. 

For the rotating braneworld with a tidal charge, we have performed the full numerical 
analysis of the epicyclic frequencies, focusing on the stable orbits, where the radial epicyclic 
frequency reaches its maximum value. Here we have found that while the large enough value 
of the positive tidal charge is not supported by the observations of high frequency QPOs, 
the situation is different for the negative tidal charge. For the large enough value of the 
negative tidal charge, the braneworld black hole exhibits close similarities with the Kerr 
black hole, including the appearance of the approximate 2 : 1 ratio between the vertical 
and radial epicyclic frequencies. We have also found that in the case under consideration, 
the over-rotating braneworld black hole and the extreme Kerr black hole are almost in- 
distinguishable in the high-frequency QPOs observations. We note that these results are 
also substantially new. On these grounds, we concluded that the potential existence of the 
over-rotating braneworld black hole with the negative tidal charge does not confront with 
modern observations of the black hole systems in the universe. 
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TABLE I. Orbital and vertical frequencies at ISCOs (M = IOMq) 





direct orbits 


retrograde orbits 


a 


/ M 


Tms/M 


Vo[Hz) 


ve[Hz) 






ve[Hz) 







00 


6 


00 


217.73 


217.73 


1.00 


6.00 


217.73 


217.73 


1 


00 





10 


5 


67 


235.32 


231.91 


1.01 


6.32 


202.54 


205.15 





99 





20 


5 


33 


255.93 


248.03 


1.03 


6.64 


189.28 


193.91 





98 





30 


4 


98 


280.49 


266.52 


1.05 


6.95 


177.59 


183.79 





97 





40 


4 


61 


310.32 


287.97 


1.08 


7.25 


167.20 


174.65 





96 





50 


4 


23 


347.48 


313.16 


1.11 


7.55 


157.92 


166.33 





95 





60 


3 


83 


395.41 


343.20 


1.15 


7.85 


149.56 


158.74 





94 





70 


3 


39 


460.42 


379.58 


1.21 


8.14 


142.02 


151.81 





94 





80 


2 


91 


556.00 


423.99 


1.31 


8.43 


135.12 


145.38 





93 





90 


2 


32 


721.41 


474.68 


1.52 


8.72 


128.84 


139.47 





92 





99 


1 


24 


1348.05 


304.86 


4.42 


8.99 


123.19 


134.10 
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TABLE II. Frequencies at direct stable orbits (r > rjsco , a = M = lOMg 



r/M 



2.00 
2.30 
2.90 
3.50 
3.80 
4.10 
4.30 
4.60 
4.90 
5.00 
5.30 
5.60 
5.90 
6.20 
6.50 
6.80 
7.00 



i/r(Hz) i^e(Hz) i^o(Hz) 



234.08 

244.60 

237.48 

217.25 

206.12 

195.14 

188.05 

177.86 

168.28 

165.23 

156.49 

148.36 

140.81 

133.80 

127.30 

121.27 

117.48 



484.35 

462.27 

398.45 

337.58 

311.02 

287.13 

272.60 

252.71 

234.89 

229.37 

213.94 

200.05 

187.51 

176.17 

165.87 

156.50 

150.71 



835.85 

712.99 

538.85 

423.96 

380.61 

344.02 

322.69 

294.50 

270.12 

262.72 

242.40 

224.53 

208.73 

194.67 

182.11 

170.83 

163.93 



T^e/J^r T^o/J^r ^^o/^d 



2.07 
1.89 
1.68 
1.55 
1.51 
1.47 
1.45 
1.42 
1.40 
1.39 
1.37 
1.35 
1.33 
1.32 
1.30 
1.29 
1.28 



3.57 
2.91 
2.27 
1.95 
1.85 
1.76 
1.72 
1.66 
1.61 
1.59 
1.55 
1.51 
1.48 
1.45 
1.43 
1.41 
1.40 



1.73 
1.54 
1.35 
1.26 
1.22 
1.20 
1.18 
1.17 
1.15 
1.15 
1.13 
1.12 
1.11 
1.11 
1.10 
1.09 
1.09 
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TABLE III. The highest radial frequency and the associated vertical and orbital frequencies at 
direct orbits (M = IOMq) 



n li\/r 
a,/ ivi 


r / A/T 
' max / 






7/ (\ir,\ 








0.00 


8.00 


70.71 


141.42 


141.42 


2.00 


2.00 


1 


00 


0.10 


7.58 


75.73 


151.17 


152.60 


2.00 


2.02 


1 


01 


0.20 


7.15 


81.52 


162.37 


165.68 


1.99 


2.03 


-I 

1 


Oz 


0.30 


6.70 


88.27 


175.45 


181.26 


1.99 


2.05 


1 


03 


0.40 


6.24 


96.29 


190.74 


199.97 


1.98 


2.08 


1 


05 


0.50 


5.76 


105.99 


209.23 


223.29 


1.97 


2.11 


1 


07 


0.60 


5.26 


118.03 


231.78 


252.90 


1.96 


2.14 


1 


09 


0.70 


4.71 


133.53 


260.59 


292.71 


1.95 


2.19 


1 


12 


0.80 


4.11 


154.57 


298.90 


349.83 


1.93 


2.26 


1 


17 


0.90 


3.42 


185.95 


353.95 


442.93 


1.90 


2.38 


1 


25 


0.99 


2.45 


243.45 


447.52 


662.11 


1.84 


2.72 


1 


48 


1.00 


2.42 


245.34 


450.40 


671.62 


1.84 


2.74 


1 


49 



TABLE IV. The highest radial frequency and the associated vertical and orbital frequencies at 
retrograde orbits (M = IOMq) 



a/M 




z/r(Hz) 


//e(Hz) 


^/o(Hz) 








0.00 


8.00 


70.71 


141.42 


141.42 


2.00 


2.00 


1.00 


0.10 


8.41 


66.31 


132.83 


131.72 


2.00 


1.99 


0.99 


0.20 


8.81 


62.42 


125.20 


123.24 


2.01 


1.97 


0.98 


0.30 


9.21 


58.95 


118.38 


115.74 


2.01 


1.96 


0.98 


0.40 


9.60 


55.84 


112.25 


109.08 


2.01 


1.95 


0.97 


0.50 


9.98 


53.03 


106.70 


103.10 


2.01 


1.94 


0.97 


0.60 


10.36 


50.48 


101.65 


97.72 


2.01 


1.94 


0.96 


0.70 


10.73 


48.16 


97.04 


92.85 


2.02 


1.93 


0.96 


0.80 


11.10 


46.03 


92.81 


88.41 


2.02 


1.92 


0.95 


0.90 


11.47 


44.07 


88.92 


84.35 


2.02 


1.91 


0.95 


1.00 


11.83 


42.26 


85.32 


80.63 


2.02 


1.91 


0.95 
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TABLE V. The highest radial frequency and the associated vertical and orbital frequencies for the 
positive tidal charge 





direct orbits: 


/3 = 0.1M^ M 


= lOMo 






n IM 








i/o(rlzj 








0.00 


7.81 


72.78 


145.69 


145.69 


2.00 


2.00 


1 


00 


0.09 


7.41 


77.83 


155.54 


156.99 


2.00 


2.02 


1 


01 


0.19 


6.99 


83.65 


166.83 


170.17 


1.99 


2.03 


1 


02 


0.28 


6.56 


90.42 


179.91 


185.77 


1.99 


2.05 


1 


03 


0.38 


6.12 


98.43 


195.30 


204.58 


1.98 


2.08 


1 


05 


0.47 


5.66 


108.10 


213.71 


227.80 


1.98 


2.11 


1 


07 


0.57 


5.17 


120.05 


236.29 


257.38 


1.97 


2.14 


1 


09 


0.66 


4.64 


135.38 


264.85 


296.79 


1.96 


2.19 


1 


12 


0.76 


4.07 


156.08 


302.77 


353.12 


1.94 


2.26 


1 


17 


0.85 


3.39 


186.74 


357.31 


444.84 


1.91 


2.38 


1 


24 


0.95 


2.42 


244.13 


452.24 


666.79 


1.85 


2.73 


1 


47 



TABLE VI. The highest radial frequency and the associated vertical and orbital frequencies for 
the negative tidal charge 





direct orbits: 


fi = - 


M'\ M 




1OM0 






a/M 


^max / 


z/,.(Hz) 




/^o(Hz) 


VQjVr 






0.00 


9.64 


56.20 


112.29 


112.29 


2 


00 


2.00 


1 


00 


0.14 


9.09 


60.69 


121.06 


122.28 


1 


99 


2.01 


1 


01 


0.28 


8.53 


65.94 


131.25 


134.11 


1 


99 


2.03 


1 


02 


0.42 


7.95 


72.15 


143.25 


148.35 


1 


99 


2.06 


1 


04 


0.57 


7.36 


79.63 


157.61 


165.82 


1 


98 


2.08 


1 


05 


0.71 


6.74 


88.84 


175.16 


187.83 


1 


97 


2.11 


1 


07 


0.85 


6.09 


100.51 


197.15 


216.53 


1 


96 


2.15 


1 


10 


0.99 


5.40 


115.87 


225.74 


255.80 


1 


95 


2.21 


1 


13 


1.13 


4.64 


137.33 


264.91 


313.73 


1 


93 


2.28 


1 


18 


1.27 


3.77 


170.50 


323.56 


411.81 


1 


90 


2.42 


1 


27 


1.41 


2.53 


236.74 


431.78 


663.00 


1 


82 


2.80 


1 


54 
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